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A strong J/ip suppression in central Au+Au collisions has been observed by the PHENIX ex- 
periment at the Relativistic Heavy Ion Collider (RHIC). We develop a hydro+ J/ip model in which 
hot quark-gluon matter is described by the full (3+l)-dimensional relativistic hydrodynamics and 
J/ip is treated as an impurity traversing through the matter. The experimental J/ip suppression 
pattern in mid-rapidity is reproduced well by the sequential melting of Xc, ^ > and J/ip in dynami- 
cally expanding fluid. The melting temperature of directly produced J/ip is well constrained by the 
participant-number dependence of the J/ip suppression and is found to be about 2T C with T c being 
the pseudo-critical temperature. 

PACS numbers: 25.75.-q, 25.75.Dw, 25.75.Nq 



I. INTRODUCTION 

A new state of matter composed of deconfined quarks 
and gluons, the Quark-Gluon Plasma (QGP), is expected 
to be formed in relativistic heavy-ion collisions if the sys- 
tem reaches a temperature larger than the critical value 
T c ~ 160 — 190 MeV as predicted by the lattice Quan- 
tum Chromodynamics (QCD) [l|. To find the experimen- 
tal evidence of QGP, the heavy quarkonia (J/ip, ip', Xc, 
and T) have long been considered as the most promising 
probe, since they are supposed to melt away due to the 
color Debye screening at sufficiently high temperature Q • 
Recent lattice QCD studies show that J/ip would survive 
even up to about 2T C while Xc and ip>' would be dissoci- 
ated at different temperature 0,0, Q. Therefore, the 
heavy quarkonia may be used as a thermometer of QGP 
in relativistic heavy- ion collision experiments [||. 

Recently, high statistics data of Au+Au collisions at 
the center of mass energy per nucleon (^/snn) of 200 
GeV at the Relativistic Heavy Ion Collier (RHIC) in 
Brookhaven National Laboratory (BNL) become avail- 
able ■ It is observed that J / ip yield in central Au+Au 
collisions at RHIC is suppressed by a factor of 4 at mid- 
rapidity and 5 at forward rapidity relative to that in p+p 
collisions scaled by the average number of binary col- 
lisions. Cold nuclear matter (CNM) effects due to the 
gluon shadowing and nuclear absorption of J/ip at the 
RHIC energy were evaluated from the J/ip measurement 
in d+Au collisions Trend of the J/ip suppression 
in d+Au collisions as a function of rapidity is reasonably 
well reproduced by the gluon shadowing with the Eskola- 
Kolhinen-Salgado (EKS) parameterization [9(1 and the 
nuclear absorption cross section cr a b s ~ 1 mb [8j. Then, 
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a comparison of the J/ip yield observed in Au+Au col- 
lisions at RHIC to that expected from CNM effects [l(| 
reveals that J/ip is anomalously suppressed in central 
collisions (Ti| . 

So far, two scenarios have been proposed for the 
J/ip suppression at RHIC energies. The first scenario 
is based on the sequential melting of the charmonia [l2l ] : 
The J/ip suppression previously observed at Super Pro- 
ton Synchrotron (SPS) in European Center for Nuclear 
Research (CERN) may be ascribed to the complete melt- 
ing of %c and ip 1 and associated absence of the feed down 
to J/ip, while at RHIC, extra melting of J/ip could be 
seen due to higher initial temperature [13|. This ap- 
proach is still in a qualitative level and the space-time 
evolution of the system has not been taken into account. 
The second scenario is based on the substantial charmo- 
nium dissociation due to gluons and comovers supple- 
mented with the regeneration of J/ip due to cc recombi- 
nation [3, 15, 10. [17| . In this scenario, a large amount of 
regeneration is expected at RHIC in comparison to SPS. 

The purpose of this paper is to make a first attempt 
to investigate the sequential charmonia suppression in a 
dynamically evolving matter produced in Au+Au colli- 
sions at RHIC energy. We develop a hydro+ J/ip model 
in which hot quark-gluon matter is described by the full 
(3+l)-dimensional relativistic hydrodynamics [l8| and 
J/ip, Xc, and ip' are treated as impurities traversing 
through the matter. We will focus on the J/ip data in 
the mid-rapidity region throughout this paper, since the 
hydrodynamical description of various observables is best 
established in that region. The J/ip suppression in for- 
ward rapidity will be remarked at the end of this paper. 

Hereafter we study the survival probability of J/ip, 
Sj°^ = Raa/R ( aa A , as a function of the number of par- 
ticipants A part [11]. Raa is the standard nuclear modi- 
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fication factor defined by 



Raa 



dN 



Au+Au 



/dy 



(N col )dN"li/dy 



(1) 



dN^t/dy 



are the invariant 



where dN^J +Au /dy and 

J r /ip yield in Au+Au and p + p collisions, respectively. 
(N co i) is the average number of nucleon-nucleon colli- 
sions. i?^A M 1S a contribution to Raa originating from 
the CNM effects constrained by the data of the d+Au 
collisions. 



II. HYDRO+J/V MODEL 

A. Full 3D Relativistic Hydrodynamics 

We solve the equations of energy-momentum conser- 
vation d^T^" — in full (3+l)-dimensional space-time 
(T,x,y,r] s ) under the assumption that the local thermal 
equilibrium is reached at initial time to(=0.6 fm/c) [lj|. 
Here r and r) s are the proper time and the space-time 
rapidity respectively, x and y are transverse coordi- 
nates. The centers of two colliding nuclei are located 
at (x, y) = {b/2, 0) and (—b/2, 0) before collision with an 
impact parameter b. The ideal hydrodynamics is charac- 
terized by the energy-momentum tensor, 



= (e + P)u>*u v - Pgf- 



(2) 



where e, P, and are energy density, pressure, and lo- 
cal four velocity, respectively. We neglect the finite net- 
baryon density which is small near the mid-rapidity at 
RHIC. The equation of state (EOS) of massless parton 
gas composed of u, d, s quarks and gluons is employed 
in the QGP phase (T > T c =170 MeV), while a hadronic 
resonance gas model is employed for T < T c . The ini- 
tial energy density distribution in the transverse plane 
is parameterized in proportion to d 2 N co \/dxdy according 
to the Glauber model. The pseudo-rapidity distribution 
of charged particle multiplicity dN ch /dri for various cen- 
tralizes observed at RHIC has been well reproduced by 
the (3+l)-dimensional hydrodynamics simulations with 
the above setups [Hj]. 

The same space-time evolution of the QGP fluid ob- 
tained as above, which is now open to public [201 ] . has 
been also exploited to study hard probes such as az- 
imuthal jet anisotropy, nuclear modification factor of 
identified hadrons, and disappearance of back-to-back jet 
correlation [2lJ. Figure Q] shows the local temperature as 
functions of x (upper) and y (lower) for various proper 
time at mid-rapidity (j] s = 0). The left and right pan- 
els are for the case at b — 2.1 fm (A part = 351) and 
b = 8.5 fm (A/ part = 114), respectively. 
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FIG. 1: Local temperature in the unit of T c as functions 
of the transverse coordinates, x (upper) and y (lower), for 
various proper time (r) at mid-rapidity (r/s = 0). The im- 
pact parameter b — 2.1 fm (the left panel) corresponds to 
(Afpart) = 351, while b = 8.5 fm (the right panel) corresponds 
to (iVp art > = 114. 



B. J/ip inside QGP 

We distribute the initial J ftp's in the transverse plane 
at rj s = according to the spatial distribution of A co i 
calculated from the Glauber model for a given impact 
parameter. Since we focus on the mid-rapidity, only 
transversal motion of J ftp is considered. Transverse mo- 
mentum (pt) of J ftp is distributed according to the in- 
variant p T spectrum of J/ip measured by PHENIX 0]. 
This corresponds to a weak coupling case where J/ip is 
assumed to follow the free-streaming path inside the hot 
matter unless it is dissociated [Hj]. 

The survival probability of J/ip suffering from disso- 
ciation along its path inside the expanding fluid may be 
expressed as 



Sj/^{xj/^{t)) = exp 



T Ais {T{ Xj/ 4T')))dT' 



, (3) 



where Tdi s {T(xj^(T))) is the decay width at tempera- 
ture T and xj/^(t) is the transverse position of J/ip at 
proper time r. Information on rdi s (T) is rather scarce at 
present because the matter is still in non-perturbative 
region T/T c ~ 1 — 2. We will first make a simplest 
threshold ansatz according to the similar one in [23| : 
r dis (T > Tj/^) = oo and r d i s (T < Tj^) = where T J/4> 
is the melting temperature of J/ip. A more sophisticated 
parametrization of Pdi s (T) will be discussed later. 
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We assume the same p T distribution for Xc and tp' as 
that for J/tp and estimate their survival probabilities in a 
similar way as Eq. ([3]) with a common melting tempera- 
T x . Thus, the total survival probability 



ture Tjj,/ = 



of J/tp is obtained by taking into account the total feed 
down fraction / FD from \ c + tp' to J/tp as 



r*tot 

J/i> 



(1 — fpn)^J/ip + /fd^Xi 



(4) 



where S x denotes the survival probability of %c or tp' . 

Although the precise values of 7j/^ Xc ,^/ should be 
eventually calculated from lattice QCD simulations with 
dynamical light quarks, we take those as free parameters 
to fit the present experimental data of Sjyy,- The feed 
down fraction / FD is also treated as a free parameter. 
At present, its value measured at different energies has a 
large uncertainty, 15% — 74% [24j. Also it may well de- 
pend on collision energy and has not yet been measured 
at the RHIC energy. As we will see below, the shape 
(magnitude) of S^fy as a function of N paYt is essentially 
determined by T j/ip (/ FD ). 



C. Numerical Results 

Experimental survival probability of J ftp in the mid- 
rapidity region in Au+Au collisions at RHIC is shown as 
filled circles in Fig. [2J where gluon shadowing and nuclear 
absorption with cr a b s = 1 mb are taken into account as 
CNM effects [ll|. Bars and brackets correspond to the 
uncorrelated and correlated errors with respect to iV part , 
respectively 0. Boxes correspond to the uncertainties 
associated with the nuclear absorption cross section of 
mb - 2 mb [H. 



best fit parameters (x 2 = 0.86 for A part > 50), 

(T JN „ T X J FD ) = (2.02T C , 1.22T C , 0.30), (5) 

while Sj/jf, and S x are shown by dashed and dotted lines, 
respectively. 

Decreasing of S^h with increasing A part is reproduced 
quite well with the scenario of sequential melting in ex- 
panding fluid: Onset of genuine J/tp suppression around 
A^part ~ 160 can be clearly seen by the dashed line in 
Fig. [2] which results from the fact that the highest tem- 
perature of the matter reaches to Tj at this centrality 
(See Fig. Q]for comparison). Gradual decrease of 
above A par t ~ 160 reflects the fact that the transverse 
area with T > Tj/^ increases. 
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ber of participants Apart , where melting temperature of J/tp is 
1.90T C (dashed), 1.96T C (dotted), 2.02T C (solid with open cir- 
cles), 2.08T C (dash-dotted), and 2.14T C (solid). 
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FIG. 2: Survival probabilities Sf/^ (solid line), Sj/$ 
(dashed line), and S x (dotted line) in the hydro+ J/tp model 
as a function of the number of participants A par t with 
{Tj /4 ,,T Xl f FD ) = (2.02Tc,1.22T c ,0.30). Filled symbols are 
the experimental survival probability in the mid-rapidity of 
J/tp in Au+Au collisions at RHIC [fl,ITl|]. 



The solid line corresponds to the net survival probabil- 
ity Sjyy, obtained from our hydro+J/-0 model with the 



Sensitivity of S^K to the melting temperature of 
J/tp is shown in Fig. [3l Here, / FD and T x are fixed 
to be 30% and 1.22T C , respectively, as given in Eq. 
The shape of the theoretical curve of 5%, is sensitive to 
Tj/^p: Even the 6% change of the melting temperature 
of J/tp from the best fit value 2.02 T c causes a substan- 
tial deviation from the experimental data. To make this 
point clearer, we show, in Fig. [U the \ 2 contour plot 
in the plane of Tj/^/T c and T x /T c with / FD being fixed 
to 0.30. The cross symbol corresponds to Eq. ([5]) which 
gives minimum \ 2 ■ Solid and dashed lines correspond to 
the lcr and 2cr contours, respectively. It is seen that Tj 
can be determined in a narrow region around 2T C , while 
T x is not well determined because the feed-down effect is 
not a dominant contribution to Sjfi^. 

Tj/,,p is also insensitive to the nuclear absorption 
cross section. Taking cr a b s = mb (2 mb) as an 
upper (lower) bound [Tl|, we found (Tj/^, T x , / PD ) = 
(2.00T C , 1.02T C ,0.35) with minimum X 2 of 0.97 for 

A par t 

minimum \ 



> 50 ((T J/4 ,,T x ,f PD ) 
2 of 0.87 for K 



(2.02T C ,1.02T C ,0.15) with 
part > 50). The change of cr abs 



affects only the overall normalization of Sfh, whose ef- 
fect is mostly absorbed by the change of / FD . In other 
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words, the value oiTj/^, is intimately linked to the shape 
of Sf^ around A part ~ 160. 

It is noticeable that the RHIC data analyzed with the 
state-of-the-art hydrodynamics leads to a rather stable 
value for the melting temperature of J/tp to be around 
T/T c ~ 2 [25|. In fact, this number is not in contra- 
diction to the result of recent quenched and full lattice 
QCD simulations which suggest the survival of J/tp above 
T c [1,0, Hi- To make a quantitative comparison, however, 
we need to wait for further progresses in lattice QCD sim- 
ulations with dynamical quarks [5|. 
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FIG. 5: Survival probability S^/^ as a function of iVp ar t, 
where melting temperatures of J ftp and x an d feed down 
fraction are 2.02T C , 1.22T C and 30%. a), b), c), and d) are 
the results using the decay width parameterized by Eq. © 
with a of 0.1, 0.2, 0.3, and 0.4 GeV, respectively. 



change from the shoulder shape to the monotonic shape 
as a increases. The data are not reproduced well for 
a > 0.1 GeV. It is unclear whether the gap between the 
data and the solid curve for a > 0.1 GeV can be filled 
by the charmonium regeneration, since the regeneration 
is a monotonically increasing function of -/Vp art . 



E. J/tp in Forward Rapidity 



D. Thermal broadening of J/tp 

The results presented so far are obtained by assum- 
ing that the decay width of charmonia changes abruptly 
from zero to infinity across the melting temperature. To 
make smooth connection with the scenario of the J/tp 
dissociation due to thermal quarks and gluons in QGP, 
we introduce the following parametrization for the char- 
monia width motivated by the T-dependence in [26| : 



r dis (T < T m ) = a(T/T c - l) 2 



(6) 



with rdis(T > 2jm) = oo. a is nothing but the thermal 
width of J/tp at T/T c = 2. The J/tp dissociation width 
in NLO perturbative QCD calculation suggests a > 0.4 
GeV [26(, while we leave it as a free parameter to take 
into account the fact that the system is still in the non- 
perturbative regime. 

Figure [5] shows the survival probability as a function of 
-'Vpart, where the model calculations shown in a), b), c), 
and d) are obtained for a — 0.1,0.2,0.3, and 0.4 GeV, 
respectively. The melting temperatures of J/tp and x 
and feed down fraction are fixed to be 2.02T C , 1.22T C and 
30%. It is seen that the suppression pattern has a smooth 



Let us here make a brief remark on the J/tp suppres- 
sion at forward rapidity which is reported to be even 
stronger than that at mid-rapidity [7[. Experimentally, 
there is still a sizable systematic error for relative nor- 
malization between the two cases [ll| . This could re- 
duce the apparent difference of Raa between the data 
at mid and forward rapidities. Theoretically, a part of 
the difference may come from gluon saturation. Accord- 
ing to Color-Glass-Condensate (CGC) model, the charm 
production at forward rapidity (y ~ 2) may become 
as low as 60% of that at mid-rapidity in most central 
Au+Au collisions [27J. If we take this number literally, 
the J/tp suppressions in both rapidities become almost 
comparable. There is also a possibility that <7 a b s depends 
on the J/tp rapidity and is larger in forward rapidity [28| . 
With all these uncertainties, we leave the J/tp suppres- 
sion in forward rapidity as an interesting subject for fu- 
ture works. 



III. SUMMARY 

In summary, we have investigated J/tp yield at mid- 
rapidity in the relativistic heavy ion collisions by using 
hydro+ J/tp model which serves as a dynamical approach 
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to the sequential charmonia suppression. The space-time 
dependence of the local temperature T of the fluid is de- 
scribed by the state-of-the-art relativistic hydrodynamic 
simulations, and the J/ip dissociation is assumed when 
T exceeds the melting temperature Tji^. Comparing 
the survival probability S^L ) from the hydro+J/^ model 
with the one in Au+Au collisions at ^/s NN = 200 GeV 
obtained by the PHENIX experiment at RHIC, observed 
J/ip suppression pattern is described quite well with 
Tj/^/T c ~ 2. This value is determined primary by the 
Afpart dependence of S^K and is insensitive to the mag- 
nitude of the nuclear absorption. Occurrence of the sup- 
pression of directly produced J/ip around 2 T c is in accor- 
dance with the spectral analysis of J/ip in lattice QCD 
simulations. Trend of the suppression can be described 
well with the small decay width inside the fluid below the 
melting temperature. 

There are a number of issues to be studied further: 



They include (i) better treatment of the J/ip propaga- 
tion in hot matter beyond the free-streaming or the com- 
plete thermalization, (ii) detailed study on the azimuthal 
anisotropy (1*2) of J/ip , (hi) effect of relative velocity be- 
tween J/ip and the hot fluid on the pt distribution of 
J/ip [29l | . (iv) better understanding of the suppression 
in forward rapidity in connection with the color glass 
condensate, (v) the application to lighter system such 
as Cu+Cu one, (vi) effect of the charmonia regenera- 
tion processes to SyJ^ in conjunction with the sequential 
suppression, and (vii) the theoretical investigation of the 
melting temperatures and the widths of J/ip, Xc and ip 1 
in full QCD simulations. 
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